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■ Abstract 

' The local reggeon field theory is studied perturbatively taking advantage of the PT symmetry 

in the Hamiltonian formulation. In the lowest non trivial order we show that the pomeron 
C^l interactions renormalize the slope. In the same order we find a non local pair potential acting 

(— I ' between pomerons, which has a singular structure. However the analysis of the scattering 

operator shows that at small coupling constant bound states do not appear so that the two- 
i particle spectrum is not changed. 

> ■ 

o : 

. 1 Introduction 
(N ■ 

■ In recent years the study of strong interactions in the so called Regge kinematics has experienced 
! a rebirth of interest due to advances in the Quantum Chromodynamics studies, originated from 

QQ I the seminal work about the BFKL [1] pomeron. Subsequent investigations have unveiled the 

■ basic effective interaction [21 |3l 3] of propagating BFKL pomerons. Equations to resum tree 
^ . diagrams at large Nc were derived for onium- nucleus scattering [5l E] (BK equation) and nucleus- 

I nucleus scattering [7]. An effective quantum field theory describing the pomeron interaction 

^ ■ at large Nc was constructed in [H E], which in principle allows to find quantum corrections 

. , related to pomeron loop diagrams. Similar studies were conducted in the alternative dipole 

and JIMWLK techniques, where also attempts to sum pomeron loops were made under certain 
drastic approximations [9]. 

Pomeron loops play a secondary role in the scattering with nuclei at not too high energies. 
But for the proton scattering and at asymptotic energies their contribution cannot be neglected. 
Summing pomeron loops is a formidable task. So, as a first step, it is worthwhile to study some 
features of a much simpler quantum field theory of the local supercritical pomerons, introduced 
by V.N.Gribov many years ago to sum reggeon diagrams which describe interacting pomerons in 
the phenomenogical approach. This local reggeon field theory (LRFT) possesses many features 
similar to the QCD pomeron theory. In particular, for the supercritical pomeron with the 
intercept a(0) > 1, it is also non-perturbative at high energies, so that loop contributions 
have to be summed by some technique, which may prove to be useful also for the QCD. Some 
very beautiful results in this direction were obtained for the even simpler LRFT living in zero 
transverse dimension ("a toy model'), which in fact reduces to the quantum mechanics. In 
particular it was shown that loops indeed play a decisive role at high energies and essentially 
transform the initially supercritical pomeron into a weakly subcritical. These results were in 
fact obtained long ago [TOl E] and were recently re-analyzed in [121 131 El • Unfortunately 
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it is not straightforward to extend these findings to the reahstic LRFT living in two transverse 
dimensions. Differentg approximationshave lead different authors to either predict a phase 
transition or claim complete inconsistency of the model. More studies are needed to clarify the 
situation. 

In this paper we draw attention to the fact that in the Hamiltonian approach LRFT belongs 
to the class of models with a non-Hermitean Hamiltonian, which possesses a certain symmetry 
equivalent to the PT symmetry in the ordinary quantum mechanics. This was noted in [14J for 
the toy model in zero transverse dimension but remains true also for the realistic case. The PT- 
symmetric non-Hermithean Hamiltonians have been extensively studied for some time |15^ I16j 
and we are going to apply some of the techniques developed in this study for the LRFT with 
a supercritical pomeron. In particular using the PT invariance we demonstrate that the non- 
Hermitean Hamiltonian of LRFT has only real eigenvalues and transform it to an Hermthean one 
by a suitable similarity transformation. As for zero transverse dimensions in jl4j . we are able to 
do this only by perturbation expansion in the coupling constant of the pomeron interaction A. 
This does not allow us to find the true asymptotic of the amplitudes in LRFT at high energies 
with all loops taken into account but gives results valid up to rapidities y ~ when terms 

up to order are included. Thus, at small A, higher orders in the perturbation expansion of 
the Hamiltonian allow to study higher and higher rapidities. 

In the course of our study we have to perform renormalization of the loop correction to the 
pomeron intercept, which is divergent in the lowest order, as is well known. As a result we find 
that the slope of the pomeron trajectory is decreased by interaction. We find that this is the 
only change in the supercritical pomeron spectrum at order A^. 

The plan of the paper is as follows. In the next section we briefly review some features of 
PT symmetric quantum systems relevant to our study. Then we introduce the LRFT model 
and its Hamiltonian formulation and discuss The PT and CPT symmetries of the Hamiltonian. 
Next two chapters are devoted to the calculation of the metric operator and the Hermithean 
Hamiltonian in the lowest non-trivial order. In the fifth section we study and renormalize the 
single particle term of the Hamiltonian and derive the effective two particle non-local potential. 
In the sixth section we study the scattering states associated with this potential and show that 
no bound states are present at perturbative level. We then draw our conclusions. Some more 
technical details have been added in a few appendices. 

2 PT symmetric Quantum Mechanics and QFT 

Since it was noted [15] that there exist non hermithean Hamiltonians H having a real spectrum 
bounded from below, provided boundary conditions for the wave functions of the associated 
Sturm-Liuville differential problem are properly defined, a lot of investigations have been done. 

An important observation has been that a class of such Hamiltonians possesses an unbroken 
PT symmetry, so that eigenstates of H can be chosen to be eigenstates of PT. After investigating 
the scalar product for the state space it has been found that there is a natural choice (/, 5) = 
/ (ix[Pr/]*(x)5((x), t denoting the transpose operation, which gives an associated norm conserved 
in time. The problem of this choice is that it is not leading to an Hilbert space with a positive 
norm. The states are splitted in two classes with positive and negative PT-norms. The good 
news has been that it is generally possible to find a new symmetry operator, denoted by C, 
whose eigenvalues are precisely the sign of the PT-norm. 

Having [C, PT] = and [C,H] = 0, one can define a new scalar product by {f\g) = 
/p dx[CPT fY'{x)g{x), which leads to a positive norm conserved in time and therefore to a physi- 
cally acceptable probabilistic interpretation. In such a case one can also define observables to be 
operators such that O* = CPTOCPT, relation which coincides with the usual hermiticity condi- 
tion in the conventional quantum mechanics where C = P. We stress that the main conceptual 
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point here has been that if one insists to define a quantum theory with conserved probabihty 
using a PT-symmetric non-Hermithean Hamiltonian one has to use the new mathematically 
well defined CPT-scalar product, which depends on the Hamiltonian itself. 

On the other hand in many physical situations the scalar product of the wave functions is 
well defined from the start and cannot be changed. With this scalar product the probabilities 
given by the norm of the wave functions will not be conserved in time. However this does not 
present any difficulty, since either they are not conserved physically, as for decaying channels, 
or evolution in fact goes not in time but in rapidity, as in LRFT, when requirement that the 
norm be conserved is absent. Thus the new norm introduced in earlier studies is not relevant 
for LRFT. However other points have a direct application. 

In the next section we shall find that the Hamiltonian of the LRFT indeed posseses the PT 
symmetry and has real eigenvalues. Construction of the new symmetry operator C depends 
on the Hamiltonian. If one knows eigenstates of H then one can use them to explicitly find a 
representation of C. A more convenient way, which admits a perturbative approach [17], is the 
following. Consider a system whose Hamiltonian has the form 

H = Ho + XHj, (1) 

where the free part is given by a Hermithean Hq and the interaction part by an anti-Hermithean 
Hj. Define the parity operator P with = 1 to transform H into H'^ , which implies 

[Ho,P]=0, {Hj,P} = 0, P^ = l. (2) 

Now one looks for the symmetry operator C satifying 

[C,H] = [C,PT]=0 (3) 

in the form 

C = e^P (4) 
where Q is an Hermithean operator. From ^ we find a relation [17J 

2Xe^Hi = [e^,H], (5) 

which can be solved perturbatively, using the expansion Q = XQi + X^Qs + by requiring 

[Ho,Q,] = -2Hj, [Ho,Q3] = -^[[Hi,Qi]Qi] (6) 

and so on. 

The operator e*^ can be used to define a similarity transformation which maps the PT- 
symmetric Hamiltonian H onto an Hermithean Hamiltonian h with the same set of eigenvalues. 
It is easy to show that also 

g-Q^gQ = (7) 

Indeed we have e^PH — He^P = Multiplying by P and using PHP = we find ([7|). As a 
result, we find 

h = e-«/2i7e«/2 = eQ/2^tg-Q/2 ^ ^g) 

so that h is the equivalent Hermithean Hamiltonian. Once Q is known as a power series in A 
the Hamiltonian h can also be found in the same form: h = + A2/i(2) + A^/i^^) + ... where 
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(9) 



and so on. 

In the next sections we show that it is possible to apply these general results to the LRFT 
and find the symmetry operator Q and the Hermithean Hamiltonian h in the first non-trivial 
order of the perturbation theory. 



3 PT symmetry and operator Q in the LRFT 

The LRFT can be defined as a theory of two fields x) and x) depending on rapidity y 
and transverse coordinates x with a Lagrangian density 



fi — a 



'v2)0 + a0t(x) 



(10) 



where /i > is the intercept minus unity and a' is the slope of the pomeron trajectory. With 
// > the corresponding functional integral is divergent and the only way to define the theory 
beyond the set of perturbative Feynman diagram is the analytic continuation from /x < when 
the theory is well defined. A constructive way to do this continuation is the Hamiltonian 
approach. One sets up a quasi-Schroedinger equation for the wave function 



dy 



-H-fiy), 



where the Hamiltonian has the form 



H = Ho + XHi 

with the free part given by 

Ho = J (fx{-^x^\x)^{x)+a'V^\x)V(l){x)), 

the interaction part 

Hi = i d^x(j)\x) (j)^ (x) + (j){x) (pix), 



and the standard commutation relations between 6 and 



{x),cp\x')]=6\x-x'). 



(11) 



(12) 



(13) 



(14) 



(15) 



The scattering amplitude with the target ('initial') state ^i{yi) at rapidity yi and the projectile 
( 'final') state ^/(?/2) at rapidity y2 > yi is defined as 



iAfi{y2 - yi) = (^/(2/2)|e-^(^^-^i)|^.(yi)). 



(16) 



One can demonstrate that the perturbation expansion in powers of A of this expression repro- 
duces the standard Reggeon diagrams of the LRFT and also that (fT6]) satisfies the requirement 
of symmetry between the target and projectile (see [4]) 

The Hamiltonian is not Hermithean and our first task is to demonstrate that its energy 
levels are all real. This is of course trivially seen in the perturbation theory. Since Hq has its 
eigenstates with a fixed number n of pomerons, the energy change can only be accomplished 
by action of an even number of interactions Hj. However we can also prove it on more general 
grounds. To this end we consider symmetry operations applied to H. In correspondence with 
the definitions in the previous section we introduce parity P as a transformation of the fields 



4>{y,x) ^ -()){- x), (l)^{y,x) ^ -(p^-x). 



(17) 
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It follows that indeed PHP = and of course P^ = 1. Next we introduce the 'time reflection' 
T as taking the complex conjugate of all coefficient functions without changing the fields. 

It is evident that P and T commute and that their product PT will leave both parts of the 
Hamiltonian intact. So we indeed find a PT symmetry of the LRFT Hamiltonian 

[PT,H] = 0. (18) 

As a result, if an eigenfunction of H is presented as the action of some operator depending on 
(j)^ on the vacuum \I'o, then it has to be of the form 

F(i0^)^o, with F{z*) = F*{z), (19) 

i.e. F has to be a real function of its complex argument. The eigenvalue is then 

(l-li^l^) = {^o\F{-i(j))HF{i(l)^)\ipo) 

and is obviously real, since the Hamiltonian can also be written as a real function of icf)^ and 
and the total mumber of operators (p and (j)^ has obviously to be the same. 
Following the technique discussed in the previous section our aim is to find an appropriate 
Q-operator, which will allow us to pass to a Hermithean Hamiltonian h by transformation ([8|). 
Once we find the latter, the amplitude will be given by 

iAf,{y2 - yi) = (e<3/2M/^(y2)|e-'^(j^^-^i)|e-«/2M/,(yi)). (20) 

Thus evolution will be accomplished by the Hermithean h and operators e^*^/^ will transform 
(differently) the initial and final states. In particular the pomeron Green function at rapidity y 
and momentum k will be given as as 

S\k - k')G{y, k) =< 0\(l){k)e^/'^e-y^e-^/^(l)^{k')\0 > . (21) 

To construct Q we shall use perturbative equations ([6]) presented in the previouss section. 

In this paper we shall restrict ourselves to the first non-trivial order in the triple pomeron 
coupling constant A, that is constructing Qi and h^'^\ To accommodate to our notations in [1] 
we seek Qi in the form 

Qi = -2-^y"d^2;i(i2x2d^X3(^/i(xi,a;2,a;3)(/)|(/)2(/'3 + /2(a;i,X2,a;3)0|0^(?i3) , (22) 

where we denote cpi = (p{xi) etc. To calculate the part of the commutator [Qi,Hq] which 
contains function /i we need to know 

[(p\4>24>3, (l)Hx)(l){x)] = 5^(X2 - x)(/)[(/)3(/)(x) + 5^(X3 - x)</)[ (/)20(x) - (5^(xi - x)(l)\x)(j)2(t>3 

and 

[0l</'203,V</.t(x)V</.(x)] 
= V6^{X2 - x)0|(/.3V0(x) + V6^{X3 - x)(t>\(t)2V(t>{x) - V6^{xi - x)V(t>^x)^2^3 , 

where V refers to differentiation in x. So this part of the commutator [Qi, Hq] takes the form 

• „ 3 

[Q^^\Ho] = -2- l[d^xi4(t>2h(- f^ + c,'i-2Vl + Vj))fi{xi,X2,X3). (23) 
^ i=i 

Here it has been taken into account that /i(a;i, X2, X3) is symmetric in X2 and X3 
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In full analogy we calculate the commutators related to the part with f2- 

[44(1)3, cl)\x)cl){x)] = 6\x3 - x)44(/.(x) - 6\x2 - x)4{x)4Mx) - 5\xi - x)0t(x)403 
and 

[014^3, V</.t(x)V</.(x)] 

= V5^{x3 - x)4i42^(l>{x) - V6'^{x2 - x)v4{x)4Mx) - V5'^{xi - x)V0^(x)4</'3 , 
which gives the second part of [Qi, Hq]: 

■ . 3 

[Qf\H^\=2- y n^Wi4'/'3(-M + a(-2V? + V^))/2(xi,a;2,X3). (24) 



i=l 



It has been taken into account that f2{xi-,X2,X3) is symmetric in xi and X2- 

To satisfy the first of the conditions ([6]) which determines the form of Qi we have to require 



and 



(-^ + q'(-2V2+V2))/i(xi,X2,X3) = - ^x5\xl - X2)5^ {xi - X3) (25) 

( - /i + a'(-2V? + V^)) /2(xi, X2, X3) = fi5^{x3 - xi)52{x3 - X2) . (26) 

These equations are trivially solved in the momentum space. We define the Fourier trans- 
forms by 

fi{ki,k2,k3)= J Y[[cfx,e-"'^^^)fi{xi,X2,X3), i = l,2. (27) 



Then we find 



and 



( - ^ + a'i2kl - kl))h{kuk2, ks) = -^i{2T,f6^{ki + k2 + k^) (28) 



(^-^ + a'{2kl - kl))f2{ki, k2, k^) = fi{27rf6\ki + ^2 + k^) . (29) 
So taking into account the symmetry properties of the functions /i and /2 one has 



{2Trf6{ki +k2 + fcg) 
'l,-a'{kl + kl-kl] 



fiiki, k2, ks) = n- ,,, 2 , ,,2 — (-^^) 



and 

rn ; 1 \ {2Trf6{ki+k2 + k3) f 1 1 \ ^qi^ 

f2{kl,k2, k3) = -/^ ^_^/(fc2^fc2_fc2) = -/l(^3, ^2-^1) • (31) 

Note that the denominators in (j30p and (jSip may vanish. The requirement that Qi be a 
Hermithean operator implies that the singularities in (j3U|) and (|31|) be circumvented from op- 
posite sides or taken both in the principle value prescription. The study of the transformed 
Hamiltonian in the next two sections reveals that at order there appears a pairwise inter- 
action of pomerons similar to the pairwise interaction of normal non-relativistic particles. In 
the latter case the interaction potential is standardly real, which leads to invariance under time 
reflection. The pomerons are not propagating in real time, but rather in the imaginary one cor- 
responding to the rapidity. Still, as mentioned above there exists a similar invariance consisting 
in changing signs of (/> and 4 s-nd taking complex conjugate of the coefficient functions. This 
requires the pair potential to be real as in the normal theory. As we shall see in the next section 
this requirement requires that the functions /i and /2 be real, so that their singularities be taken 
in the principal value sense. This circumstance will be implicitly understood in the following. 
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4 The transformed Hermithean Hamiltonian h^'^^ 

In this section we shall find the second order Hamiltonian /i^^^ determined by second of Eqs. 
dH). We present the interaction term Hj as an integral over three momenta qi, i = 1,2,3 using 

^{x) = J ^e^'^m (32) 

and similarly for (ft^x). With this normalization the fields will obey the standard commutation 
relations in the momentum space 

[(t>{k),(l)^k')]=S^{k-k'). (33) 
The interaction Hamiltonian acquires the form 

i=l 

+ {2TT)H\qi +q2- q3)cl>HQi)^HQ2mq3)] • (34) 
In the same manner we get operator Qi as 

Ql = -2- / n ^2, h)cl>Hk,)<l){k2mks) + /2(-fcl, -k2, h)cl>\ki)cl>\k2)<P{h)) . 

^ i=l ^ 

(35) 

To find the second order Hamiltonian h^^^ we have to calculate the following 4 commutators: 

C2 = [cPHk^)cl>{k2)cl>{ks),cp\q,)cl>Hq2)HQ3)] , 
Cs = [^Hh)(t>Hk2)Hk3)AHqi)HQ2)Hq3)] , 

C4 = [<t>Hk,)4>Hk2)Hk3),4>Hqi)<PHq2)4>{q3)] ■ 

They are all trivially found: 

Ci = -5\q2 - k^)<pHqi)<P{q-3)<P{k2)<P{k3) - S^qs " A:i)</.t(<7i)</'(<72)0(A:2)</'(A:3) 

+ s^iqi - k2)^Hki)Hk3)4>{q2)Hq3) + 5\qi - k3)(l)Hki)4>{k2)4>{q2)Hq3) , (36) 

C2 = -5\ki - q3)(PHqi)(PHq2)Hk2)(p{k3) 
+6{k2 - q2)^Hqi)^Hkl)Hk3)Hq3) + S{ks - q2)(l)\qi)^Hkimk2)4>{q3) 
+ S{k2 - qi)^Hkimk3)(l)\q2mq3) + 5(^3 " qi)^Hkimk2W {q2)Hq3) , (37) 

C3 = -C2{ki,k2,h ^ Q'l, 92,93) 
= 5\qi - h)^Hkl)(l)\k2)ct>{q2)Hq3) " S{q2 - k2)(l)\ki)4>\qi)<l,{q^)cl){h;) 
-5(93 - k2)<i>Hki)<p\qi)<P{q2)^{k3) - <5(92 - k^)<p\qMq3)<p\k2)<i>{k^) 

- 5(93 - k^)cpHqi)^{q2)<p\k2)<P{k^) , (38) 
C4 = Cl{ki,k2,k^ ^ 93,92,91) 
= -6\k2 - q3)^Hqi)(p\q2)(pHki)Hk3) - S\ki - 93)0^ (91)0^ (92)<^^(A:2)0(A:3) 
+ S\k3 - q2)(t>Hki)4>Hk2)<PHqi)Hq3) + ^'(^3 - 9i)0^(^i)</'^(^2)</'^(92)</'(A;3) . (39) 
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Passing to /i^^^ we get the following 4 terms 

1 f (Pk2(fk^d?q2 



(27r)2 J IX- a'{kl + kl- {k2 + ks^) 
(f{k2 + k3)ct>{ks)cl>{q2)4>{k2 - q2) - <t>\q2 + ^2 + k^)<t>{q2)(t>{k2)<t>{k^)) , (40) 

,(2) 1 [ A d'^^jd^li (27r)"^"(gi + 92 - g3)(27r)^(5(A;i - ^2 - fcg) ^ „^ , , , . 
=--2] n T^^^p ^-a-(fci + fci-fc^) ^^1^^' 



1 1 /• d-k2<fk:',(Pqi 



■{24>\qi)4>\k2 + k3)cp{ks)cf){qi + ^2) 



2 (27r)2 y ;u - a'(A;| + A;| - {k2 + A;3)2) 
+ 2<^t(^2 + k3)cl>{k3)(l>H(li)H(li + ^2) - 4>\qi)4>\k2 + ks- qi)cf>{k2)Hk3)) , (41) 

(2) ^1 f A 6?2fc,d2g. (27r)2(52(gi - 92 - g3)(27r)2(5(fci + ^2 - fcg) ^ „^ , , , . 
^3 =^2j}l^^ ^-a'{kl + kl-ki) C,{k„k2M<lu<l2,<ls) 

.-a^^^^^+W ^ ^^^^^^^^^^^^ ^ 

+ 2(/.t(fc2 + q3)<i>{q3)<t^\ki)(t>{ki + ^2) - 0^(A;i0^(fe2)</.(A;i + A^s - 93)</'(g3)) , (42) 



,(2) ^1 [ A ^'fe^^'g^ (27r)'(52(gi + 92 - g3)(27r)^(5(fci + ^2 - fcs) ^ 

^4 =+2yn-(^;^)2 ^-a-(fc2 + fc2_fe2) C,{k,MM\quq2,q.) 

^ l_ f (fk2(fk:id^q2 

(27r)2 I Pi- a'{kl + kl- {k2 + k^Y) 

(f{k2 - 92)0^(^2)0^(^3)0(^2 + fea) - 't>\ki)ct>\k2)(t>Hq2)'i>{k2 + fe3 + 92)) . (43) 
We have obviously 

so that the total second order Hamiltonian h is Hermithean. 

5 Single- and two-particles terms in h and renormalization 

All terms in /i^^) gpji^ jn^o the pomeron number conserving, hf^ +hf \ and pomeron number 

(2) (2) 

changing: h\ , with AA'" = —2 and /14 with AN = +2. The two latter terms will contribute to 
energy levels only in the order and can be neglected in the order with we restrict ourselves 
here. 

The pomeron number conserving terms may be rewritten in the normal form. As a result 
we get two contributions with two or four field operators, describing single- and double- particle 
parts. The single particle part is found to be 

7(2) 2 ^ /■ (fk2(fk3 ,4.,, , N ,/, , N 

W = -J2^2^^ J ^,-a'{ki + kl-{k2 + ksr/ ^''' + ^^^^(^^ + 
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= J (fk(l)^{k)(j){k)A^^h{k), (44) 

where 

is the shift in the pomeron energy in order A^, so that the total pomeron energy is 

e{k) = -fi + a'k'^ + X^A^^h{k). (46) 
One easily finds that A^'^^e{k) can be presented (see Appendix A) as 

A(2)e(fc) = +^ r ^e-(/^WfeV2) . (47) 
47ra' Jo X 

The integral obviously exists only for < —a'k'^/2 and diverges at x = 0. To regularize it we 
require that at A: = the change of pomeron energy vanishes, that is e(0) = — In a way this is 
a definition of the renormalized intercept. One can obtain the same result choosing the standard 
dimensional regularization approach. With this condition the regularized Ae(A;) is 

A('hreg{k) = ^L.(,+a'kV2) _ 

^ 4:iTa' Jo X ^ ' 

1 , / olk^\ , 
■In 1 + ^— . (48) 



47ra' V 2/i 

It exists for any vahics of /i and is real for /i > 0. Taking into account that the initial pomeron 
trajectory is determined only up to terms linear in A;^, we approximate the energy shift as 

A(2)e,es(A;) = -^fe^ (49) 

so that the net effect of the single-particle term in is to renormalize the slope: 

1 



ol a'^^^ = a' -\^- — . (50) 



The two-pomeron interaction term h^p}^^. can be presented in the form corresponding to 

transition ki, k2 ^ qi, q2 

h% = j c^'A;id2A;2dVdW(9i + g2-A;i-A;2)F(')(gi,g2|fci,A;2)<^t(gi)<^t(g2)<^(A;i)<^(fc2), (51) 



where 



V^^\qi,q2\kiM)=-^ 



27r2 - a'{kl + {k2 - qiY - ql) 27r2 /x - a'{ql + {q2 - h)^ - k^) 

^8^,1- a'{kl +kl- (fci + k2Y) n-a'{ql + ql- {qi + gs)^) ^^^^ 

(symmetrization in qi,q2 and ki,k2 is implied). As mentioned, poles in the two terms have to 
be understood in the principal value sense for the potential to be not only Hermithean but also 
real. 

This potential is non-local and degenerate. To more clearly see its properties consider a 
case when gi + ^2 = + ^2 = corresponding to the two-pomerom exchange for the forward 
scattering amplitude. Then denoting qi = —q2 = q and ki = —k2 = k we have a potential 

V^'\q\k) 11 11 



27r2 II - a'{k'^ + {k + qf - q^) 27r2 - a'{q^ + {k + qf - k^) 
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The last two terms depend only on the initial or only on the final momenta. Since the poten- 
tial falls rather slowly at high momenta, its integration over q or k meets with a logarithmic 
divergence. (See the form of the kernel in the coordinate space in Appendix 2.) 

The effect of this degenerate pair potential is not quite clear in the general case. Considera- 
tions in the next section, restricted to the case qi + q2 = ki + k2 = 0, tell that the spectrum of 
the pomeron states will not be changed by this interaction. So in its presence the two-pomeron 
states will continue to have their total energy 

E2iki,k2)=€iki)+eik2), (54) 

with e(fe) given by ([16|) but the wave functions will become changed by the standard scattering 
operator. To the second order in A 

^k.,kM,Q2) = Sym{5\k^ - q^)5\k2 - 92) + A % , ^^^,'''^t\'''\. , , J ' (^5) 
1^ e{qi) + e{q2) - e{ki) - e{k2) ± tO J 

where symbol Sym means symmetrization in qi and (72 and signs of iO correspond to in- or 
-outgoing waves. 

For the asymptotic of the Green function we shall find from the two-pomeron states 

5\k - k') [ d2A;id2A;2(0|(/>(fc')e«/2|M/fc^,,,)e-^^^(*^^''^^)(*fc,,fc,|e-«/2^t(fc)|0) , (56) 



where the matrix elements can be easily computed by perturbations in A. This asymptotics will 
be true at large y until y ~ l/X^. Of course one also will have similar contributions from states 
with the number of pomerons greater than two, but the corresponding matrix elements will be 
of the higher order in A. 

6 The Schroedinger equation with a degenerate potential 
6.1 Problem 

Consider the Schroedinger equation in the 2-dimensional momentum space 

{e{q) - E)^{q) = - f S kV {q\k)^{k) . (57) 



To simplify notation we rescale E to exclude all terms independent of q in e{q) and have in our 
case e{q) = 2a'q^. Our potential has a structure ([53|) : 

V{q,k)=v{q)+v{k)+Vi{q,k), (58) 

where Vi{q\k) has the normal properties and vanishes as any of the arguments go to infinity. 
Our aim is to study the spectrum E of the solutions to Eq. (|57|) 

Obviously Eq. ([57j) can have both solutions corresponding to the scattering states and 
to bound states. In the former case we standardly convert this equation into the Lippman- 
Schwinger equation presenting 

where T[q\l) (T from now is no more the "time reflection" operator) satisfies 
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If this equation can be solved it corresponds to the scattering state with the incident momentum 
I and energy E = e(/) which belongs to the continuous positive spectrum. 
For the bound state iPe{q) with energy ii^ < we analogously present 

M,) - (01) 

with an equation for Ie 

Our aim is to study possible solutions of Eqs. (|6U|) and (j62|) with a degenerate potential (|58|) . 
6.2 Continuous spectrum 

Putting (j58p into (j60p and suppressing the fixed argument / in T we have 



Tiq)=v{q)+c+Vi{q,l)+v{q) d^k 



.2, Tik) 



e{l) - e{k) 



. v{k)T{k) , (.V^{q,k)T{k) 



where we have denoted the part of V independent of q 

v{l) = c. (64) 

We also denote 

.2. nk) _ f,2,v{k)T{k) 



We find an equation 

(n^ = r P ^ (^ ^ rJ'\-y,(n\ -U T/, ( n 1\ -i- f rPh- 

e(/) - e{k) 



T(q) =c + e+{l + d)v{q) + V,{q, I) + j d'k^^^^^. (66) 
Correspondingly we present 



Tiq)=c + e + {l + d)viq)+Tiiq). (67) 

The equation for Ti{q) is 

Tiiq) = ic + e)xiiq) + il+d)x2iq) + ViiqJ)+ j d^k^^^^^^, (68) 



with 



.2,. V,(q.k) 



and 



For the two constants d and e we obtain equations following from their definition 

d = (1 + d)h + (c + e)/o + / d^k-^^^^ (71) 
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and 



where 



e = (c + e)Ii + (1 + d)l2 + J (fk 



v{k)Ti{k) 
e{l)-e{ky 



d'k- 



v'^ik) 



(72) 



(73) 



e{l)-e{ky 

To solve Eq. (j68p with additional conditions ()7ip and (j72p we first solve this equation for 
three different inhomogeneous terms 



T^\q)=4\q)+ j dPk 



,,, Vl{q,k)Ti'\k) 
e(l)-e{k) ' 



where i = 1, 2, 3 and 



(74) 
(75) 
(76) 

Now we put this solution into the equations ([7T|), ([72]) to obtain a system of two linear 
equations for d and e: 



ri'''^('Z)=Xl,2(<z), 4'\q)=V^{q,l). 
Prom these three solutions we obtain the solution to Eq. (j68p as 

Ti(g) = (c + e)T['\q) + (1 + d)rf ^(g) + (g) 



d{h + J2 - 1) + e(/o + Ji) + c(/o + Ji) + A + + = , 

d(/2 + K2) + e(/i + i^i - 1) + c(/i + Ki)+h + K2 + K^ = Q, 
where /„ are defined by ([73]) and 



'e{l)-e{ky 



eil) - e{k) ■ 



(77) 



(78) 



All quantities in fact depend on the fixed momentum I. Solution of the linear system ([77p is of 
course trivial. The determinant is 



D = {h + .h- l)(/i +K1-I)- (Jo + Ji){h + K2 



(79) 



and so 



D 



(c(/i + Ki) + I^ + K2 + Ks) (Jo + /i) - (c(/o + Ji) + Ii + J2 + J3) (/i + i^i - 1)] , (80) 

(c(/o + Jl) +h + J2 + J3) (/2 + i^2) - (c(/l + Kl) + /2 + i^2 + Ks) {h + J2- 1)] . (81) 

Of course the resulting d and e are functions of the fixed momentum /. With thus determined 
d{l) and e{l) Eq. ([76p gives the final solution to the Lippmann-Schwinger problem. The only 
difficulty may generally arise in case D{1) = which may only happen at some specific values of 
I and leads to certain singularities of the scattering matrix at these values of momentum, which 
we consider improbable. 

For future reference, here we present orders in A for different quantities defined in the previous 
derivation. Obviously 

D = l + 0{X% /„~A2"; xi(2)~A2W 



J2 ~ ETi, ~ A^ i^2 ~ A^ 



(82) 
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Then it follows from ()80p and (|8ip that d ~ and e ~ and in the lowest approximation 
(order A^) the scattering matrix T is given just by the total potential V, as expected. 

However in our case there is a new problem. The constant Iq is in fact divergent. In the 
limit Iq — > oo we find that the determinant grows linearly with Iq: 

D = -Io{h + K2). (83) 

The denominators of (j80p and (|8ip also grow linearly with Iq. So in the limit /q — > oo we find 
finite values for both d and e: 

l + d = -^— ^, e + c = 0. (84) 

With these values we find in this limiting case 

Ti (q) = rf ) (q) - I^T'f ^ (<?) • (85) 

So again the solution in all probability exists but the constant e is automatically adjusted to 
exclude the constant term c = v{l) from the original Lippmann-Schwinger equation ()66p . 

As a result, we find that even for divergent Iq the Lippmann-Schwinger equation has a 
solution for any positive energy, so that the spectrum is continuous and covers all positive 
values of energy 

Note that in the limit Iq ^ oo orders of d and e in powers of A are radically changed. Now 
the determinant D ~ A^ and as a result 1 + d ~ 1/ A^ and e ~ A^ . As a result already in the 
lowest order A^ the scattering matrix Ti acquires additional terms from T^^^: 

Ti{q) = V^iq\l)-^X2{q). (86) 

Turning to the full scattering matrix (j67p we find that it acquires a term of the order unity 

T(,|0 = -^+r.(,). (87) 

It may be considered as a renormalization term for the scattering matrix in the limit /q ^ oo. 
Of course appearance of this term is due to the implicitly made assumption that A^/q >> 1 as 
/q — > oo. Different relations between the small A and large Iq will lead to different results. 

To illustrate the described procedure for the solution of Lippmann-Schwinger equation in 
Appendix C we calculate the scattering matrices T and Ti up to order A^ for the pair pomeron 
potential V^^\ Eg. ([53]) . for the forward case gi + ^2 = + ^2 = 0. The found expressions are 
long and not very interesting but they show that the procedure is quite feasible and does not 
encounter any new complications. 

6.3 Discrete spectrum 

Putting (f58l) into (f62]) we now obtain 



tE{q) = V{q) I d^k^^^ + J d^k^-^ +jd^k ^ _ ^^^^ 



As before we denote 



,2, tE{k) f j2jV{k)tE{k) 



d^k ^^' e= d'k ,\ / . (89) 

E-e{ky J E-e{k) ^ ' 



14 



We find an equation 

tEi,) = e + dvi,)+Jd^k^^^}^. (90) 

We present 

tE{q) = e + dv{q)+tiE{q) (91) 

to find an equation for tiE 

tlEiq) = exiiq) + dx2{q) + J ^'^^^^^^ (92) 
wliere similarly to the continuous spectrum case 



and 



with the two conditions to determine d and e 



d = dh + elo + / cPk^^^^ (95) 
J hj — e(k) 



E-e{k) 



and 



,2,v{k)tiE{k) 



e = eli + dh + / d'k (96) 

h — e[K) 



The integrals In are the same as in (I73p with e{l) — > E. 

Obviously the solution to Eq. (j92p can be presented as a sum 



tii5(g) = et!2((?) + dtfi((?), (97) 



(1 2) 

where the two functions tj^^ satisfy 



4fa) = X>(<,)+/Al%M|pl, , = 1,2 (98) 

After functions are known, one finds a homogeneous system of linear equations to 

determine d and e: 

d = elo + dh + eJi + dJ2, e = eh + dh + eKi + dK2, (99) 
where now, similarly to ()78p . 

All the coefficients in the system (|99p depend on i?. The value of the bound energy E is found 
from the condition of existence of solutions to the system (j99p : 

(II + J2- l)(/l + i^l - 1) - (/o + Jl){h + i^2) = 0. (101) 

Now consider the case /q ^ oo. Then Eq. (jlOip reduces to 

h{E) + K2iE) = 0, (102) 

which determines a possible bound state in this limiting case. 

As we have seen, at small values of A h is of order A"^ and K2 is of order A^ Since h{E) 
cannot vanish (it is strictly negative for < 0) Eq. (ll02p cannot be satisfied. So for small values 
of A there are no bound states when Jq ^ 00. 
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7 Conclusions 

We have generalized the technique of constructing an Hermithean Hamiltonian for the PT sym- 
metric LRFT model, developed in [4J for the toy model in zero transverse dimensions, to the 
realistic case of two transverse dimensions. The complexity of the latter model makes both 
the derivation and analysis of the found Hamiltonian not straightforward already in the lowest 
non-trivial order in the coupling constant A of the triple pomeron interaction. In particular the 
divergence of the pomeron intercept has to be eliminated by renormalization. Also the found 
pair interaction between pomerons is both singular and degenerate. It has required a separate 
study of the Schroedinger equation with degenerate potentials, which may have a wider scope 
of applicability. 

As a result we have found that at small A the total impact of the pomeron interaction at 
order A^ is reduced to the change of slope. These results allow to study the asymptotic of any 
scattering amplitude at large rapidities y in the region 



However in the course of our study it became clear that at finite A the pomerons may form 
bound states, whose presence will drastically change this asymptotics. 

Note that the experimental values for a' and A are roughly a' ~ 0.25 (GeV/c)~^ and A ~ 0.33 
(GeV/c)^^, so that the actual parameter of the perturbative expansion is A^/a' ~ 0.4, which is 
not so small. So to estimate the possibility to apply our results to the realistic processes one 
has to study the NNLO (terms of order A'^). If they happen to be relatively small then one can 
study possible pomeron bound states using our potential and, say, the variational methods. 

We stress that the perturbative study of our Hermithean Hamiltonian cannot give the true 
asymptotic of the theory at y — > oo, which remains unperturbative. To find it one has to search 
for non-perturbative techniques to construct this Hamiltonian. 
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A Single pomeron energy shift 

We have to calculate the integral 



We present the J-function as an integral over r and the denominator at < —a'k'^/2 as 



1 « y « {a'/X^f 



I = 




(103) 




(104) 



to get 




(105) 



Integrals over k2 and k-^ give the same result: 




(106) 
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So we find 

I = A r dx I d\-\.e-"''^e^^^^^'^"-^^ . (107) 
A Jo J {xa'Y ^ ' 

Next we integrate over r 



y d^ri^'-^ = 27rxa'e-^"''='/2 (iqs) 



to finally find an integral over x 



1 = -^ r ^e-(^+"''=V2) . (109) 

B The coordinate space potential between pomerons 

The two-pomeron potential is non-local in the coordinate space. We shall limit ourselves with 
the case qi + q2 = ki+ k2 when the in the momentum space the potential is given by (j53p . Then 
in the coordinate space its kernel is defined as 

V^^\y\x) = f f^e-'^y+"'''V^^\q\k) . (110) 



. (2vr)4 

The non-degenerate part of the kernel comes from the first two terms in (j53p and is 
We rewrite the first term in this expression as 



27r2 J (2^)4 ' ' ^ ' /X - a'(A:2 + ^2 _ ^2) 

= - , \ , f (fq(fk(fKd^re-'^y+'^''+'''^'^+^-^^ \ , ^ . (112) 

{2-Kfa' J ^ m + ak^ + hq^ + CK? ^ ' 

where a = c = — 1,6 = 1 and m = n/a' We further present 

1 



and do the Gaussian integrals in k, q and k assuming that they exist, that is for positive a, b 
and c. Transition to the desired values of a, b and c will be achieved by analytic continuation. 
We get 

Vl^\y\x) = --4^^ r ^e-^^ [ d2^g-(:.+r)V4a5-(y-r)V46€-r-V4c€ _ (1^4) 



{2Tr)^a' abc Jo 
The exponent in the integrand in r has the form 

-iPr^ -2sr + xya + yyb)/i^, 

where 

P = l/a + l/b + l/c, s = {-x/a + y/b) . 
so that after the integration over r we find 
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where 

.2 Jl „,2 



Integration over ^ gives 

where zi = 2\/mi Inserting the desired values a = — 1, 6 = 1 and c = — 1 we find abcP = — 1 so 
that 

z^ = J'-Hl^ (118) 
V a' 

and thus 

4'^(yk) = ^^^Ki(.i) (119) 

If zi is real then the potential is falls exponentially. If zi is pure imaginary, then putting 
z = —i\z\ and taking the real part according to the principal value prescription the potential is 
expressed via the oscillating Neumann function: 

^/?(2/l-) = I^^^N,(|zi|) (120) 

The second part of the potential Vi is obtained as 

vS\y\x) = V}?{-X\ -y) = ^^1ki(z2) , (121) 

where 



,, = J'J^^^^. (122) 



The degenerate part of the potential is simpler. We have 



We calculate the first part by taking jj, = —2a'v and assuming u > to subsequently continue 
to negative u. So 

The integral in k is trivially done to give 

Continuing to negative u and putting V vx"^ = —i\/\v\x^ we finally find 



Taking into account the requirement of hermiticity, the second part of the degenerate potential 
will be given by 
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C Lowest order scattering matrix for the pomeron interaction 

As in Section 4 we restrict ourself to the forward case 91+^2 = ^1 + ^2 where the full potential 
V{q, k) is given by ([53]) . so that Vi{q\k) is given by the first two terms in ([53|) and 

Obviously in our case Iq is divergent. We assume that it is regularized in some manner (say 
restricing integration by A; < A). As A ^ 00 /q — > 00 logarithmically. This leads us to the 
expressions (fSHI) and (f87l) for Ti and T respectively valid up to order A^. To calculate these 
scattering matrices all we need is to find X2iQ) and I2, keeping in mind the reality of the 
potential and the prescription given in Eq. (I60p . 
Calculation of I2 is of course trivial. We find 

A"^ 1 r 1 /, N 1 1 , ^ 

/2 / = -rr— 7-^72 2 72 2 0^-?^-^^+—' 

where we use a convenient notation 

= (130) 

Calculation of X2{q\1) is a bit more complicated. It consists of two terms coming from the 
two terms in the potential Vi. The first term can be written as 

(1) A"^ r cfk 



USn^a'^ J [{k + - p^]{k^ - m?W ' l^) 

= ^ j A{f)- A{m^)] , (131) 



where 



and 



^^^^^^ ! [{k + q/2Y-p^W-P)- ^^^^^ 
This latter integral is conveniently calculated using the Feynman parametrization to finally give 

^ ' ^/A {P-m^-q^/2 + VA){P-m-^-VA) ^ ' 

where 

A = (/2-m2)2-g2;2 (135) 
and in (jl34p it is assumed that A > 0. For A < we find 

Re^ r) = ^= arctg , ^ ' arctg , - . (136 

^ ^ ^ ^Za ^ V^A J ^ ^ 



7r2 



ImA(/^) = --=0(A). (137) 



The imaginary part is 



The expression for A{m^) is simpler due to cancellations between p'^ and m^. It is real: 

A(m'^) = -—arctg— . (138) 
qm 2m 



19 



The second term in x can be written in a manner similar to (jl3ip : 

y(2)f.|n = [ = ^ (B(f\ B{m^)\ 

^^'^ rnvrV^y (qk + g2_,„2)(^2_^2)(p_^2) 1287r4a'3 ^2 - ^2 V"^^' ^ ^V'^ ) ) . 

(139) 

where now ^ 

^^^'^"/(kq + g2_^2)(p_^2)- (140) 
Integration over the azimuthal angle gives 

where 

Ai = (g2 _ ,n2)2 - g2fc2 (;L42) 

and km is defined by the condition that Ai > 0. Here we have used the principal value pre- 
scription which tells that at A < the real part of the azimuthal integral is zero. Subsequent 
integration over k gives 

B(i^)^Re-^ln ^+|^;-'";| -i^O(A.). (143) 

VA2 ^^2-^2-771^1 VA2 

where 

A2 = {q^ - rr?f - q^f (144) 
and A2 > 0. For A2 < i?(/2) is real and given by 

Calculation of B{m?) obviously gives the same expressions (|143p and (jl45p in which /2 ig 
substituted by m?. This finishes calculation of X2('i'IO- 
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